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Recent developments have made possible the computa- 
tion of equilibrium dynamical correlators in quantum impu- 
rity problems. In many situations however, one is rather inter- 
ested in correlators subject to a non equilibrium initial prepa- 
ration; this is the case for instance for the occupation prob- 
ability P{t) in the double well problem of dissipative quan- 
tum mechanics (DQM). We show in this paper how to handle 
this situation in the framework of integrable quantum field 
theories by introducing "boundary interactions changing op- 
erators" . We determine the properties of these operators by 
using an axiomatic approach similar in spirit to what is done 
for form-factors. This allows us to obtain new exact results for 
P{i); for instance, we find that that at large times (or small 
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Strongly correlated electron systems, which have been 
the subject of intense study in recent years, pose many 
theoretical challenges due to their essentially non pcrtur- 
bative nature. The greatest progress have been realized 
for systems that are one dimensional, in particular quan- 
tum impurity problems. Fixed points and their vicinity 
can then be investigated with techniques of conformal 
field theory (CFT) Q, and the behaviour between fixed 
points with the Bethe ansatz g], g. Even with these 
powerful tools, the understanding of dynamical proper- 
ties is still quite incomplete. For instance the Bethe 
ansatz provides thermodynamic properties, and, with 
some effort, DC properties, including out of equilibrium 
M; but time and space dependent correlations, though 
accessible in principle [p[, have remained largely unde- 
termined until recently. It is highly desirable to make 
progress in that area, in particular in view of recent ex- 
citing experiments |g], 

In the last two years, building on formal works on in- 
tegrable massive quantum field theories (QFTs) ^, [ p^ , 
[ pT[ , it has become possible to determine some of these 
correlators exactly (ie with arbitrary accuracy, all the 
way from short to large distances) |12]. The results ob- 
tained sometimes exhibit striking behaviour, illustrating 
the difficulty to build an intuition for non perturbative 
systems. For instance, in the double well problem of dis- 
sipative quantum mechanics, it was found that the tran- 
sition from coherent to incoherent regime, if based on 
the two spin correlation function, takes place at a value 
of the dissipation a ~ \ |12| , |ia| , in contrast with the 
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value Oi = J expected before ]1J], |15|]. This result later 
raised questions about which quantity should be used to 
describe this transition M . 



Unfortunately, correlators that can be computed in the 
framework of QFT are not always the physically relevant 
ones. This is because natural objects in QFT are vac- 
uum vacuum Green functions, while in many experimen- 
tal settings, what is measurable are time evolutions for a 
system prepared in a given state that is, in fact, not even 
an eigenstate. A well known example of this is the oc- 
cupation probability P{t) in the double well problem of 
DQM |]1^. In that problem, the system has a non equi- 
librium initial preparation: for negative times, the spin 
is held fixed in the state Sz — ^, with an equilibrated 
environment. At time t — Q, the constraint is released, 
and the dynamics starts from P(0) = 1 with a factorized 
spin-environment initial state. The behaviour of P{t) has 
been the subject of much work lately Q, ||l8[. Several 
possibilities for measuring this quantity experimentally 
have also been proposed [|9|, [g^l- On the other hand, 
no exact result for P{t) has been available up to now, 
except at a = I [ fi4[ . 

In this letter, we introduce a method to compute corre- 
lators in prepared initial states that are not eigenstates. 
The key ingredient is a generalization of the bound- 
ary conditions changing operators of CFTs |l|, pl| to 
theories with boundary interactions. We illustrate the 
method by the computation of P{t). The final results 
are too bulky to be presented here; they have the follow- 
ing interesting features. For <? > ^, the behaviour of P{t) 
is incoherent, and P{t) has the general form 
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where ^Un — 1 , and Tb is an energy scale characterizing 
the tunnel splitting of the free system (see below). The 
transition from coherent to incoherent regime takes place 
a,t g = i and is related with the transition from repulsive 
to attractive regime in the sine-Gordon model. For g < 
i, P{t) has a form similar to Im) but the sum involves 
also (known) complex arguments in the exponential. At 
large times, the dominant behaviour is 
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which in turn leads to the prediction for the ratio of the 
damping factor to the period of oscillations 
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All these results are new, except for the fact that the 
transition occurs right at 5 = ^, which was demonstrated 
using perturbation around 17 = i in |17t . 

To start, we recall the concept of boundary conditions 
changing operator, and we introduce their generalization, 
in the simple case of the Ising model (closely related to 
the double well problem at g — ^). We consider the 
massive Ising model with an inhomogeneous boundary 
magnetic field 
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where apF is the usual massive free Majorana fermion 
action, a is a boundary fermion satisfying a^ = 1, as 
is the boundary spin operator, which coincides with 
^ ("0 + V') fl- If the mass vanishes, the theory is confor- 
mal invariant in the bulk. Conformal invariant boundary 
conditions are then obtained with h = (free spins) or 
h = ±00 (fixed spins). As shown by Cardy ||2^], the situ- 
ation where boundary conditions are conformal invariant 
but inhomogeneous can be represented by the insertion 
of conformal operators on the boundary. For instance, 
the case with h = ior y < and ft, = 00 for y > is de- 
scribed by the insertion of the operator $12 of dimension 
A = ^ at y = 0. 

Let us now consider a situation where conformal in- 
variance is broken: first suppose that the bulk mass is 
non zero, and that the boundary magnetic field is uni- 
form and equal to ha- This can be handled using the 
integrability of the problem Q, Q. The hamiltonian 
(for time evolution along the y direction) is diagonal- 
ized using a basis of multiparticle states (here, simply 
fermions) that have factorized scattering both in the 
bulk and at the boundary. We use in the following ra- 
pidity variables parametrizing energy and momentum as 
e — mcosh0,p — msinhO; the asymptotic states have 
then the general form \9i^ . . . ,6n)a- S and R matrix ele- 
ments relate these states with others, where some of the 
rapidities have been switched, or had their sign changed. 
Using the known S and R matrices ||9|, correlators can 
then be computed ||l2[ , ^^ , provided one knows the ma- 
trix elements of the operators (form-factors) in the mul- 
tiparticle basis iQ . For instance, for the energy density 
in the Ising model, the only non vanishing form-factors 
are 



(eiw.w)) = ^ /'gi?,(0)cosh0e-2™-inhe^ (5) 



and 
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where ti acts on functions of many variables by chang- 
ing the sign of 0„ and Rbi9) = i tanh (f - f ) ^-^^150' 

The case where the theory is still massless in the bulk is 
the most interesting for applications jl2] ; it can simply be 
obtained as a massless limit of the previous description, 
with 771 ^ and 9 —>■ ±cx). 

Now suppose that the bulk mass is non zero, and in 
addition the boundary magnetic field is inhomogeneous: 
h ^ ha for y < and h = hi, for y > O.We now have two 
different hamiltonians to diagonalize, giving rise to two 
sorts of multiparticle eigenstates. Since these multiparti- 
cle eigenstates provide complete sets of states, the scalar 
products 
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have to be non zero in general, even for disjoint sets of 
rapidities: they are nothing but the matrix elements of 
"boundary interactions changing operators" . 

Knowledge of the scalar products (0) is what is re- 
quired to compute a quantity like P(t), as we will demon- 
strate later. For the moment, let us discuss how to deter- 
mine these scalar products. Consider the simplest case, 
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Since for the Ising model, the integrable particles are 
just fermions with S — —1, G must satisfy the exchange 
relations 

G(9,,92) = -G{92,9,) = r;{9i)G{-9u92) 

= Rl{92)G{9u^92). (9) 

In addition, G must have a kinematical pole at 62 — 
9i — in. To see this, observe that there are two possible 
expressions for the one point function of the energy: for 
y > one has 
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while for y < 0, one has 
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Replace now the form factors of e by their explicit 
forms: for ( pl| ) to be the analytical continuation of ([iQ), 
we need the residue condition 



Res G{0, - in) 



-i 1 



Rai9) \ 

Rb{0)J 



fl2) 



k=l,2 



This is enough to determine the function G. We find, 
using the parametrization Ka/b = —cosh.0a/bj 
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Inspection shows that G has the same poles as those of 
Rl- Crossing in this problem reads b(O|0i,02)a = b{9i — 
ITT, 02 — *7''|0)a- There is also a non trivial monodromy: 
G{9i + 2iTT,02) = -R.aiOi)R;{9i)Gi9u02). 

Some comments are in order: the general solution for 
an arbitrary number of particles is not difficult to find 
and will be given in |23] . The previous expression depend 
only on /i^ /{,, and in fact, are only valid for hahb > 0. The 
case of fields with opposite signs requires an additional 
particle at imaginary rapidity, see |23| . 

It is interesting now to get back to the conformal case: 
let m ^ and, for instance, ha = oo and hb = 0. Setting 
e = zszp = e^ , we find from the foregoing formulas 
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where + designates up spins, F free spins. This is noth- 
ing but the form-factor of the A = rr^ operator, in agree- 
ment with |gl| . We are then led to the following spectral 
representation for the energy one point function 
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Explicit evaluation gives 



{e{w,w))- 




z\l x^ + ; 



1 



(18) 



in agreement with the result obtained in 124] in the con- 
text of conformal field theory. 

The scalar products (R) can be determined using the 
same ideas for more general integrable interacting the- 
ories. The expressions become quite bulky and won't 
be reproduced here; some very simple features determine 
most of Pit) properties, as we now explain. 

The double well problem can be mapped on a single 
channel Kondo model when the dissipation is of ohmic 
type |14]. The hamiltonian reads 
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The spins S± are spin 1/2 operators, and the Sz value 
corresponds to the two states of the dissipative model. 
The dimensionless number g characterizes the dissipa- 
tion. The quantity P{t) was defined in the introduction; 
it reads 



P{t) = (0]5,(i)]0) 
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where ]0) is the product 10)a=o ® ]+}, |0)a=o the ground 
state of the theory with A = 0. The Heisenberg opera- 
tor 5*2 (t) in ( pC| ) evolves with H\. To compute P(i), we 
insert a complete set of eigenstates of Hx on each side 
of Sz- These are multiparticle states resulting from the 
massless limit of the sine-Gordon model ]^; the spec- 
trum is richer than in the Ising case and consists of 
solitons/anti-solitons, and, for g < 1/2, of breathers 
which we denote by e = ±,1,2,... respectively. The 
determination of P{t) then requires the scalar products 
^""■■'^^(/Jp, ...,/3i]0), which are of the same nature as those 
discussed previously, together with the form factors of the 
spin operator in the massless limit. 

A subtle point has to be emphasized here: the in- 
tegrable picture is an infrared one, where the spin is 
screened; that is, it doesn't appear in the multiparti- 
cle description. As demonstrated in Q its properties 
can nevertheless be computed by using correlators of the 
current operator Q^.^ as follows 
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That the spin is up in the initial state is then taken into 
account by giving a unit charge to the state ]0); the non 
vanishing scalar products are thus those for which ^ e^ = 
1. 

At g — 1/2 the sine-Gordon model decouples into two 
massless Ising models with a boundary magnetic field 
h (X X and the previous considerations on the Ising model 
can be used. There are no bound states at that point and 
the refiection matrix for the solitons and anti-solitons is 
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Tb is describing the boundary scale and is related to A, 
the precise relation to be found in Eq. (5.7) of ]g^. The 
current operator has a simple form factor at that point 
given by (we turn to real time here) 

x{Q\d.,4>\Pi,p2)T" ^ 5,,+,,eie^^'+P-^'^ (23) 

[l + R+-{Pi)R+-{l32)\e~''^'''+^''\ (24) 

We won't need the detailed knowledge of the scalar prod- 
ucts (0); only their properties under crossing, together 



with the value of the residues at the kinematical poles, 
matter. 

Many processes contribute to the evaluation of P{t). 
The simplest takes the form 

^^ o{0mtt{f3i\0)o +(/32|9.0(t)|/?i)+. (25) 

This process gets convoluted with many others, that is 
one has to add to the factor (0|/32)J tif^'^l^) ^^^ ^^™ 
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Suppose we move the rapidity integrals to Im /3i = —in 
and Im P2 = in. Forgetting the singularities encoun- 
tered in doing so, one obtains, using crossing, the scalar 
product xif^'^lMt ~ 27r(5(/3i — /32). But one easily cheks 
that the form factor 1^ {P\dx4>\P)x vanishes, so the whole 
series adds up to sero! All what matters therefore are 
the contributions of the poles encountered in moving the 
contours. Moreover, one can show that the kinematical 
poles do not contribute; the one point function is entirely 
determined by the poles of the the reflection matrix at 
Pi 2 = log(TJ,) =pi7r/2. and one gets the well known result 
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The main feature of this computation - the fact that 
P{t) is entirely determined by the poles of the R matri- 
ces - generalizes to arbitrary values of g, since it follows 
entirely from the g-independent general properties of the 
scalar products. As a result, we can immmediately ob- 
tain interesting properties. For g > 1/2 there are only 
solitons and anti-solitons in the spectrum, for which the 
reflection matrix is independent of g, and still given by 
(P^). It follows that the behaviour is entirely incoherent, 
and that P{t) has the form given in (|l|). All terms now 
contribute since the form factor of dx(t> is non zero for 
any even number of particles at g > 2 • 

When g < 1/2 there are also m < [l/g] breathers in 
the spectrum (m integer) , with a reflection matrix given 
by[p5l 



Rm{Q) — - 



tanh( 



13-Pb 



"4(1-9) 



tanh( 



;3-/3i, 



+ i 



■King 



(28) 



with Pb — logTfc. This matrix has poles where e^ has 
non vanishing real and imaginary part: they give rise 
to oscillatory contributions to P{t), indicating that the 
behaviour is coherent in that domain. An expression sim- 
ilar to (13) can be written: the leading behaviour at large 
times (or small g) follows from the one breather, as given 
in (ph, leading to the ratio (H). This result is consis- 
tent with the expansion in. g = 1/2 — t done in pTJ. 
It also agrees with the 5 — > limit in which A ~ TrgTb/2 



p6[ . It would be very interesting to test this numerically. 
Formula (m) should be especially useful in the quantum 
optics context, where the values of g are usually quite 
small. 
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